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Gamma Ray Bursts (GRB) are the extremely energetic transient events, visible from 
the most distant parts of the Universe. They are most likely powered by accretion on 
the hyper-Eddington rates that proceeds onto a newly born stellar mass black hole. 
This central engine gives rise to the most powerful, high Lorentz factor jets that are 
responsible for energetic gamma ray emission. We investigate the accretion flow evolution 
in GRB central engine, using the 2D MHD simulations in General Relativity. We compute 
the structure and evolution of the extremely hot and dense torus accreting onto the 
fast spinning black hole, which launches the magnetized jets. We calculate the chemical 
structure of the disk and account for neutrino cooling. Our preliminary runs apply to the 
short GRB case (remnant torus accreted after NS-NS or NS-BH merger). We estimate 
the neutrino luminosity of such an event for chosen disk and central BH mass. 

Keywords: black hole physics; gamma ray bursts; MHD 

PACS numbers: 11.25.Hf, 123.1K 



1. Model of the hyperaccreting disk 

The model computations are based on the 2 dimensional, general relativistic MHD 
code HARM2D, described in Ref. 1, The nuclear equation of state is discussed in 
detail in Ref. 2, The goal of our preliminary test calculations is to investigate the 
overall structure of an accretion disk in which nuclear reactions take place and the 
gas looses energy via neutrino cooling. 
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1.1. Chemical composition 

We assume the gas to be in beta equilibrium, so that the ratio of proton to neutron 
satisfies the balance between forward and backward nuclear reactions. We assume 
neutrino cooling via electron, muon and tau neutrinos in the plasma opaque to their 
absorption and scattering. Neutrinos are formed in the following processes: URCA 
process, electron positron pair annihilation, nucleon - nucleon bremsstrahlung, plas- 
mon decay. 

1.2. Dynamical model 

We solve GRMHD equations using numerical code HARM2D. In our model, the 
gas is composed of mixture of p, n,i^e,e~ and e+, so the total fiuid energy is 
P = Po -\- u = "^niTrii -\- Ui^ where po is the total rest mass energy density and u 
is the total internal energy density and where i = p, n, i^e, e~, e+. The total gas 
pressure is P = ^ P^. All particles have a common temperature: T {Pi = UikT) and 
we assume poly tropic equation of state P = (7 — l)ii, with adiabatic index 7 = 4/3. 
We do not include the radiation and neutrino pressure in the dynamical calculation, 
but we do account for a neutrino cooling. After each timestep, we compute the 
new gas composition and neutrino cooling rates as described in Ref. 2 , The total 
internal energy is reduced using explicit method with n-sub-cycles in each time step. 
The cooling function g^, also accounts for neutrino absorption and scattering on free 
neutrons and protons. The optical thicknesses for absorptions and scatterings are 
calculated assuming approximate disk vertical thickness H ^ O.br and local proton 
and neutron number densities. Therefore in the dense disk regions the cooling is 
effectively zero. 

1.3. Initial conditions 

We start the numerical calculation from the equilibrium model based on Ref. [3] 
around a spinning black hole. The central black hole mass is SMq and the total 
mass of the surrounding gas is about 0.1 Mq. HARM2D code is designed to solve the 
magnetohydrodynamic equations in the stationary metric. We assume dimensionless 
spin of a black hole to be a* = 0.98. For higher disk masses, a fully general relativistic 
model including the evolution of the space-time and spin of the black hole should 
be considered. 

The initial magnetic field is assumed with toroidal component of vector po- 
tential ~ p (magnetic filed lines follow the constant density surfaces i.e. 
they form loops concentrated around the pressure maximum) and normalized with 

P — Ptot,gas/ Pmag 50. 

2. Results 

The torus around the spinning black hole at hyper-Eddington rates is cooled by 
neutrinos (e.g., as shown in 1-D model in Ref. Here we also study such tori. 
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however to quantify the effect of neutrino coohng in MHD simulations, we run a 
starting test model with no cooling assumed. In Figure [U we show the structure of 
such a torus (Model 1). The snapshots present the baryon density (po = riprrip + 
nnTnn^nHe^He)i g^s temperature T and magnetic {3 parameter. We also show the 
topology of magnetic field lines. 

A turbulent and dense rotating torus forms in the equatorial plane and extends 
towards the marginally stable orbit of the black hole. The polar regions are quickly 
evacuated and an electromagnetically dominated jet is formed above the poles of the 
black hole, as shown before in Refs. 5 6, The mass accretion rate is M = O.4M0/S. 

In Figure [2] we show Model 2, in which the neutrino cooling is accounted for 
in the dynamical evolution. The initial conditions in Model 2 are the same as in 
Model 1. The mass accretion rate remains the same as in Model 1, but the structure 
of the disk changes. The disk is geometrically thinner and strongly magnetized in 
comparison to Model 1. 

We compute the neutrino emissivity at t=0.02 s. The total neutrino luminosity 
corrected for absorption and scattering processes, is = 3 x 10^^ erg s~^ (Fig- 
ure [21 left panel). This number is comparable to numbers obtained from relativistic 
hydrodynamical simulations, e.g. in Refs. 0^ Also, the relativistic MHD simula- 
tions in Ref. Unreported the neutrino luminosity on the order of Ly < 10^^ erg s~^, 
depending on black hole spin (a < 0.9) and torus mass. To compute the electro- 
magnetic luminosity of the observed GRBs, one needs to consider the efficiency of 
neutrino-antineutrino annihilation process, as well as swallowing of some fraction 
of neutrinos by the black hole due to the curvature effects. 




Fig. 1. Model 1: Structure of an accretion disk in model with no cooling functions taken into 
account in the dynamical evolution. The maps show: (i) density, (ii) temperature of the plasma, 
and (iii) ratio of gas to magnetic pressure, with field lines topology (from left to right). 
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In Figure [2] (most right panel), we show the effective neutrino coohng. The 
most neutrinos are formed within lORg. We do not see a clear neutrinosphere (a 
surface at which the optical thickness for neutrino absorption and scattering is 1.0). 
Its presence as well as the effect of neutrino trapping will be accounted when the 
neutrino pressure is incorporated in the dynamical calculations. 
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Fig. 2. Model 2: Structure of accretion disk in model with neutrino cooling taken into account 
in the dynamical evolution. The maps show: (i) density, (ii) temperature of the plasma, (iii) ratio 
of gas to magnetic pressure, with field lines topology, and (iv) the effective neutrino cooling Qi, 
(from left to right). Parameters: black hole mass M = SMq, spin a = 0.98, initial magnetic field 
normalization /3 = 50, and initial disk mass M^isk = O.IMq. The snapshot is at t=0.02 s since the 
formation of the black hole. 



